Non-linear resistivity and heat dissipation in monolayer graphene 
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We have experimentally studied the nonlinear nature of electrical conduction in monolayer 
graphene devices on silica substrates. This nonlinearity manifests itself as a nonmonotonic depen- 
dence of the differential resistance on applied DC voltage bias across the sample. At temperatures 
below ~70K, the differential resistance exhibits a peak near zero bias that can be attributed to 
self-heating of the charge carriers. We show that the shape of this peak arises from a combination 
of different energy dissipation mechanisms of the carriers. The energy dissipation at higher carrier 
temperatures depends critically on the length of the sample. For samples longer than lO^m. the 
heat loss is shown to be determined by optical phonons at the silica-graphene interface. 

PACS numbers: 72.80.Vp, 72.20.Ht, 72.15.Lh, 73.40.-c 



Graphene is a promising material for electronic ap- 
plications because of its high electrical and thermal 
conductivitiesi. A limiting intrinsic room-temperature 
mobility of graphene field-effect transistors of ^2 x 
10^ cm^/Vs has been estimated^ and approached 
experimentallj«^. Its high thermal conductivitj^, far ex- 
ceeding that of copper, holds great promise for heat man- 
agement in electronics^. Fully realizing this potential re- 
quires an understanding of the scattering mechanisms in 
graphene, particularly electron-phonon scattering, which 
becomes increasingly important with increasing temper- 
ature and device currents. 

Several different techniques have been used to study 
the coupling between charge carriers and phonons in 
graphene. Measurements of the temperature dependence 
of the resistance have shown that at low temperatures, 
T< 150 K, electron-phonon scattering is dominated by 
longitudinal acoustic phonons^i^^— , resulting in a linear 
temperature dependence of the resistivity. At higher 
temperatures the resistivity becomes strongly tempera- 
ture dependent, possibly due to ripples in the graphene 
membrane^ or coupling with remote interfacial phonons 
(RIPs) at the Si02 surface^i^. This latter scattering 
mechanism has been put forward as an important heat 
dissipation mechanism in carbon nanotube and graphene 
devicesiSrJ^. Measurements of the current saturation in 
graphene lend support to the hypothesis that these RIPs 
are an important scattering mechanisnii^, although an- 
other study identifies scattering with optical phonons in- 
trinsic to the graphene lattice^^ as the dominant satura- 
tion mechanism. 

In this paper, we report electrical measurements of 
monolayer graphene field-effect transistors in the high- 
current regime at temperatures between 4 and 200 K. 
We show that overheating of the charge carriers results 
in a nonmonotonic dependence of the differential resis- 
tance on source-drain bias. By correlating this depen- 
dence with the measured temperature dependence of the 
differential resistance, we test different models of ther- 
mal transport and electron-phonon coupling, and demon- 
strate that energy dissipation is dominated by coupling 
to the remote interfacial phonons in the silica substrate. 

Graphene flakes were mechanically exfoliated onto con- 



ductive silicon substrates covered by a 300 nm thick 
Si02 layer, allowing carrier density control by applica- 
tion of a voltage between flake and substrate. Flakes 
were identified as single layers by optical contrast mea- 
surements and Raman spectroscopy, and were electrically 
contacted with Cr/Au (5/50 nm) to allow four-terminal 
measurement of the differential resistance. This resis- 
tance, dV/dl, was measured by a low-frequency (^10 
Hz) lock-in technique using sufficiently small AC current 
to prevent unintentional self-heating (typically less than 
5 nA, corresponding to an AC source-drain voltage ~100 
/iV) . The differential resistance was measured as a func- 
tion of DC voltage, Vsd, applied between source and drain 
contacts. We consider below the results of two represen- 
tative samples of different lengths, which demonstrate 
different thermal transport regimes. Sample A is shown 
in the inset of Fig. [Ub), and is 3.7 fiin wide and 16 ^m 
between voltage probes, with a typical mobility of 5700 
cm^/Vs. Sample B is a shorter sample, 4 /im wide and 
10 fim in length, with a lower mobility of 2500 cm^/Vs. 




-100 100 100 
V,(mV) T(K) 



FIG. 1: (color online), (a) Differential resistance as a function 
of DC source-drain bias measured at different bath tempera- 
tures: Tb = 4.2, 17, 26, 35, 83, and 107 K from bottom to 
top. (b) Dependence of the zero-bias differential resistance 
on the bath temperature. All measurements are taken at the 
Dirac point. The solid red line is a guide to the eye. Inset: 
optical micrograph of sample A, with the flake indicated by a 
dotted outline; the scale bar is 10 ^m. 
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Figure [Ija) shows the differential resistance of sam- 
ple A as a function of DC source-drain bias for different 
bath temperatures. These measurements were taken in 
an atmosphere of dry helium, with the gate voltage, Vg, 
tuned to the Dirac point. At the highest temperature 
shown, Tb = 107 K, a small increase in R is seen with 
increasing Vsd- As the bath temperature is decreased, we 
observe two features: the differential resistance becomes 
strongly dependent on Vsd, and a peak in the differential 
resistance is developed at zero bias. At low-temperatures 
the differential resistance behaves in a similar way to 
the temperature dependence of the zero-bias resistance 
i?o(T), which is shown for the Dirac point in Fig.[TJb). At 
low temperatures, Ro{T) exhibits a logarithmic behavior, 
due to the well known weak-localization and electron- 
electron interaction contribution a^^i^^ . The size of the 
low-temperature peak in R{Vsd) is ^ Ikfi, which is sim- 
ilar to the size of the low-temperature peak in Rq{T). 
Furthermore, the peak in R(Vsd) occurs only for bath 
temperatures below the minimum in i?(T), which oc- 
curs at Tb ~ 40 K in Fig. [IJb): it is barely resolved in 
Fig.[TJa) at Tb = 35 K, and is completely absent at Tb = 
83 K. These observations suggest that the nonmonotonic 
differential resistance originates from a hot-carrier effect. 

Further support for this hypothesis is provided by an- 
alyzing the evolution of the peak with carrier density. 
Figure ^a.) shows the change in differential resistance 
R{Vsd)~R{0) against applied source-drain bias for differ- 
ent carrier densities. (A small asymmetric component of 
the differential resistance resulting from the DC source- 
drain bias modulating the gate voltage has been removed 
by symmetrizationii.) It can be seen that the peak in 
RiVsd) is most pronounced at low carrier densities. This 
density-dependence of the peak size is also seen in the 
temperature dependence of the resistance below 50 K, 
Fig.[5Jb), suggesting this is also a manifestation oi Ro{T) 
due to a hot-carrier effect in R{Vsd)- The features in Figs. 
[T]and[2]have been seen in all monolayer graphene samples 
that we have measured. 

To make a quantitative check of the hot-carrier hy- 
pothesis, one needs to analyze the production and distri- 
bution of Joule heat in the sample. In a micron-sized 
conductor this is determined by several characteristic 
lengths: the elastic impurity scattering length I, electron- 
electron collision length l^-e, and electron-phonon relax- 
ation length le-ph- At low temperatures, inelastic pro- 
cesses are suppressed by Pauli blocking, and the relax- 
ation rates form a hierarchy, I <^ l^-e <C le-ph {I ^ 10 
nm). Estimating the electron-electron collision time in 
sample A as Te_e HEp/ikBTy ^ 0.2ns at Vq = 3V, 
one finds l^-e ~ VDr^.e Ifim. The electron-phonon 
relaxation length, which will be discussed later, can ex- 
ceed 10 fxin. In short conductors, L <ti h-e, electrons are 
scattered elastically, and no Joule heat is produced inside 
the conductor. When the conductor size exceeds several 
microns, so that Ze-e < L < le-ph, the energy of electrons 
is thermalized by electron-electron collisions, and Joule 
heat is transferred only to the electrons. This heat is 
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FIG. 2: (color online), (a) Differential resistance as a function 
of DC source-drain bias measured at Tb = 4.2K, for different 
gate voltages, Va = 0, -6, -12V with respect to the Dirac 
point (bottom to top), (b) The temperature dependence of 
the zero- bias resistance, for the gate voltages shown in (a): Va 
= 0, -6, -12V, bottom to top. In (a) the zero-bias resistance 
has been subtracted from the curves to aid comparison, whist 
in (b) Ro{T = 4.2 K) has been subtracted, (c) Differential 
resistance as a function of source-drain bias, measured with 
different lengths between source and drain contacts, taken at 
4.2 K. Measurements were taken for sample A in (a) and (b), 
and sample B in (c). 



then dissipated via diffusion of electrons into the metal 
contacts, and the resulting electronic temperature distri- 
bution is described by a universal temperature profilei^: 

T!i-)-n + i,ieVsd?{l~'-^)- (1) 

Here L is the sample length, and x is the spatial co- 
ordinate running along the length of the sample from 
—L/2 to L/2. In longer samples, L > le-ph, electron and 
lattice temperatures are equilibrated via phonon emis- 
sion. Note that in graphene the phonon contribution to 
the specific heat, Cph oc T^, dominates over the elec- 
tron contribution, Ce oc T, at all relevant temperatures. 
Therefore, nearly all Joule heat is transferred to the lat- 
tice by electron-phonon relaxation. 

To test the Joule heating picture, we first analyze the 
data for the shorter 10 fim sample. Its zero-bias resis- 
tance is shown in Fig. |31[a) as a function of bath temper- 
ature. Neglecting electron-phonon relaxation we use the 
temperature profile in Eq. ([T]) to compute the differential 
resistance: 



R{Vsd) = 2 J R{Te{x))dx. (2) 

In Fig. [3l[b)-(d) this model is compared with the differ- 
ential resistance at three different gate voltages, Vq = -3, 
-6, -12 V, with respect to the Dirac point. It can be seen 
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that Eqs. ([T]) and ([2|) describe well the central peak in 
the differential resistance, which confirms our hypothe- 
sis about its origin and demonstrates the inefficiency of 
electron-phonon relaxation in the 10 ^m sample. In tak- 
ing L in Eqs. ([ij and ([2]) to be the distance between 
the voltage probes, rather than the distance between the 
source and drain contacts, we assume that carriers are 
held at the bath temperature at each metallic contact. 
In order to test this assumption, we measured 
between the same voltage contacts in sample B, whilst 
using different contacts as source and drain such that 
the distance between the source and drain was varied be- 
tween 25 and 50 ^m, Fig.[2jc). No effect is seen in R{Vsd) 
when changing the source-drain distance, justifying our 
exclusion of any temperature difference over the interface 
between the metallic contacts and the graphene flake. 



is the electrical resistivity, and {dE/dt)ph is the energy 
dissipation rate due to phonons. The electronic ther- 
mal conductivity is given by the Wiedemann-Pranz law, 
K = TT^k^Te/3e^p. We assume that the phonon tem- 
perature is equal to the bath temperature, and impose 
the boundary conditions Te{x ~ ±L/2) = Tb, which as- 
sumes complete cooling at the leads. Given the small 
non-linearity seen in the R{Vsd) at Tb = 107 K (Fig.p, 
one might assume that up to Vsd ^ 150 mV the tem- 
perature of electrons does not significantly exceed 100 
K. Transport measurements^ have shown that optical 
phonons, either in the Si02 substrate or in the graphene 
lattice, do not contribute significantly to transport be- 
low ~200 K. Thus, we begin by considering only acoustic 
phonons in the graphene lattice. The energy dissipation 
due to acoustical phonons was found i n^^i^° : 
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FIG. 3: (a) Temperature dependence of the normalized zero- 
bias resistance of sample B at Vg= -3, -10, -16 V with respect 
to the Dirac point, (b)-(d) differential resistance as a function 
of source-drain bias at different gate voltages, Vg = -3, -10, 
-16 V, respectively. The bath temperature Tb = 4 K in (b) 
and (c), and 2 K in (d). The symbols are measured values, 
and the red lines are fits using Eq. ([l}. 



To probe electron-phonon relaxation, we have studied 
the non-linear differential resistance of a longer sample. 
Fig. H The fit based on Eqs. ^ and © (dashed line) 
is clearly inadequate, as it results in a dependence that 
is steeper than that actually observed. This means that 
Eq. ([T]) overestimates the electron temperature. The dis- 
agreement can be accounted for if we include cooling via 
phonon emission. Thus, phonon cooling is indeed im- 
portant for the 16 /im sample. We use a heat balance 
equation: 



dQ _ 

dx pL"^ 



(3) 



where Q = K{T)dTe/dx is the heat flux, k is the elec- 
tronic thermal conductivity, V"^ / pL"^ is the Joule heat, p 



dE\ 

~dr , 



DlkBE% 



(4) 



where Da = 18-21 eV is the deformation potential, 
g = 7.6 X 10^^ kgm~^ and Vs = 2x 10'' ms""' are the areal 
mass density and speed of sound in monolayer graphene^. 
The ratio of this function to the electron energy den- 
sity Ue oc Tg gives the characteristic rate of tempera- 
ture equilibration of the electrons via phonon scattering, 
which is plotted in Fig. |31[d). The nonmonotonic be- 
havior of the relaxation rate occurs due to cooling by 
acoustic phonons becoming ineffective at high temper- 
atures: since the energy of acoustic phonons is limited 
to the Bloch-Gruneisen temperature, fesIsG, multiple 
scattering events are required to equilibrate hot electrons 
{Te > Tbg)- In Fig.lljd), where = 3 V, this corre- 
sponds to Tbg = 25 K, and we see that the maximum 
in {dE /dt)ac/Ue occurs at Te ^ Tbg/'^- Also shown in 
the figure is the inverse diffusion time, r^^ = -n^D/L^^ 
which relates to the typical time for hot carriers to diffuse 
out from the sample center into the metallic contacts. In 
the 16/im sample this mechanism plays a minor role rel- 
ative to the cooling by electron-phonon scattering at all 
temperatures. 

The temperature profile in the sample can be deter- 
mined through Eqs. ((3)) and (|4]). This is plotted in 
Fig. m^e) for an applied bias of 70 mV. Using this T{x) 
for each Vsd in Eq. ([5]), the resulting R[Vsd) are plotted 
in Fig.|3{a)-(c), with Vg = 3, -6, and -12 V, respectively. 
We take Da = 18 eV for = 3 V, and = 21 eV 
for Vg = -6 and -12 V—. In comparison to purely elec- 
tronic heat transfer, the low-bias region is well described 
by the inclusion of acoustic phonon scattering, correctly 
identifying the point in Vsd where the up-turn in resis- 
tance occurs. However, at large biases, where cooling by 
acoustic phonons is inefficient, our model predicts a steep 
increase in R[Vsd), in striking contrast to the data. This 
indicates that an additional cooling mechanism becomes 
dominant at higher electron temperatures, T > 70 K. 

A plausible mechanism for carriers to dissipate heat is 
via the RIPs in the Si02 substrate, which have an unusu- 
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FIG. 4: (color online). (a)-(c) Measurements of R{Vsd) 
(black) in sample A for different gate voltages: Vg = 3, - 
6, -12 V, respectively. The smooth curves are calculations 
for the case of no phonons (dashed green), acoustic phonons 
(dotted red), and both acoustic and optical phonons (blue). 
Note that in (c) the dotted and solid curves overlap, (d) The 
rate of energy loss of electrons due to scattering by acoustic 
phonons (dotted red line), and scattering by both acoustic and 
remote optical phonons (solid blue line). The dashed green 
line shows the inverse diffusion time, (e) The calculated tem- 
perature distribution along the length of the sample for the 
three different models (same color coding) with an applied 
source-drain bias of 70 mV. The bath temperature is 4.2 K, 
and Vg = 3 V. 



ally low frequency, hiOs = 59 meV. Scattering by these 
optical phonon modes has been proposed to explain the 
temperature dependence of the graphene resistance at 
T >200 Ki^ These polar optical phonons in the substrate 
couple with carriers in the sample via a fringing electric 
field, which decays exponentially away from the substrate 
surface. The scattering by these phonons can be signif- 



icant in graphene samples due to the close proximity of 
the carriers in graphene to the underlying substrate. We 
have found the corresponding contribution to the cooling 
power—: 

where /3 = 0.025 is a dimensionless coupling constant^iil, 
which is related to the dielectric function e(w) of the sub- 
strate. Comparing this to Eq. one can see that 
optical phonons dominate the energy dissipation, see 
Fig.m^d). We have used Eq. (jSj with the two dissipation 
terms Eqs. (|H) and ([5|) to determine the temperature dis- 
tribution in the high-current regime. A steep exponential 
increase in dissipation stabilizes the electron temperature 
and this gives rise to a flat temperature profile. Fig. ^e) . 
The resulting R{Vsd) (solid blue lines in Fig. |ll[a)-(c)) 
now give a reasonable fit to the measured values over the 
entire range of source drain bias. 

To summarize, we have measured the differential resis- 
tance of monolayer graphene in response to an applied 
DC source-drain bias at low temperatures. We observe 
a peak in the resistance around zero bias, with an up- 
turn at higher biases. We demonstrate that this structure 
arises from self-heating of charge carriers. By modeling 
the behavior of the resistance as a function of bias, we 
demonstrate that surface phonons of the underlying sub- 
strate are the dominant heat dissipation mechanism for 
carriers at carrier temperatures above ^^70 K. 
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